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How large can the branching ratio of B s — > t + t be 
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Motivated by the large like-sign dimuon charge asymmetry observed recently, whose explanation 
would require an enhanced decay rate of B s — > t + t~ , we explore how large a branching ratio of this 
decay mode is allowed by the present constraints. We use bounds from the lifetimes of Bd and B s , 
constraints from the branching ratios of related b — > st + t" modes, as well as measurements of the 
mass difference, width difference and CP-violating phase in the B s -B s system. Using an effective 
field theory approach, we show that a branching ratio as high as 15% may be allowed while being 
consistent with the above constraints. The measurement of this decay, therefore, may be within the 
reach of current experiments, and can point toward a specific class of new physics models. We also 
explore the possible enhancement of this decay in models with leptoquarks and Z' , and find that in 
the latter case the branching ratio may be as much as 5%, which can alleviate the dimuon anomaly 
to some extent. 

PACS numbers: 13.20.He, ll.30.Er, 12.60.Cn 



I. INTRODUCTION 

In 2010, the D0 Collaboration reported an anoma- 
lously large CP-violating like-sign dimuon charge asym- 
metry in the B system [1, 2], which was strengthened 
by the updated measurement [3]. The weighted average 
with the older CDF results [4] gives 



A b - 



-(74.1 ± 19.3) x 10~ 4 



(1) 



which is a 3.8<r deviation from the Standard Model (SM) 
prediction A b sl SM = -(2.3 ±0.4) x 10~ 4 [5]. 

The measured dimuon charge asymmetry is a linear 
combination of the semileptonic asymmetries ag L and a SL 
in the Bd and B s sectors, respectively. Therefore, the 
new physics (NP) may contribute through either of these 
two sectors. However, since 



Ar 



q tan<^ L (q = d,s) 



where Ar g is the width difference in the B q -B q system 
and (j)^ is the CP-violating phase in the semileptonic 
decays, an enhancement in ag L is necessarily accompa- 
nied by an enhancement in tan</>^ L and/or in Ar 9 . In 
the Bd sector, since tan</>| L sa 0.075 in the SM already, a 
further large enhancement would amount to fine-tuning. 
Also, an enhancement in Ar^ would imply the NP con- 
tribution of a few percent to the branching ratios of de- 
cay modes common to Bd and Bd, which is ruled out 
by the measurements of such modes [6]. In the B s sec- 
tor on the other hand, w 0.004 in the SM [5], and 
the branching ratios of some of the decay modes, notably 
of B s — > t + t~ , have not yet been strongly constrained. 
Therefore NP that contributes to B s — > t + t~ would be a 
prime candidate to account for the dilcpton anomaly [7] . 
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Indeed, it has been shown that the only effective four- 
Fermi operators that can account for such an enhanced 
A b SL are (&rs)(tIY) and (brs)(cTc) [8]. We, therefore, 
investigate how large the branching ratio B(B S —> t + t~ ) 
can be. 

Since the desired NP is expected to contribute to B s 
and Bd decays to different extents, it may lead to a dif- 
ference in the lifetimes of Bd and B s , which is otherwise 
expected to be < 1% in the SM [5]. The recent LHCb 
measurements [9] yield T B jT Bd = 1.002 ± 0.014 ± 0.012. 
This indicates that at the 2<r level, the branching ratio 
B(B S —> t + t~) up to 3.5% is still allowed, even if there 
is no NP contribution to Bd decays. If NP contributes 
to Bd decays, this bound will be further relaxed. Note 
that this would be a large enhancement: the value of this 
branching ratio in the SM is ps 7 x 10~ 7 [10]. 

Measurements from other modes of the form b — > 



ST^T 



like Bd 



Bd 



Kt + t and Bd 



(2) K*t + t , can restrict the NP contribution to B s 



since they involve the same effective four-Fermi opera- 
tor. However the experimental information on such de- 
cay modes is very poor. The only direct bound available 
at present is the 90% BaBar limit of [11] 



B(B^ 



")l,»>14.23GeV» <0.33%. (3) 



It was suggested in [10] that a dedicated experimen- 
tal analysis of the LEP data using the missing en- 
ergy spectrum as in [12] could be used to constrain 
B(B d -> X s t+t-) and B(B S -> t+t-) to 5%; how- 
ever such an analysis has not been performed to our 
knowledge. The bound on B(Bd X s t + t~) obtained 
in [13] also weakens considerably when the theoretical 
uncertainties in B(B —> Kiv + anything) are taken 
into account. The charm counting in Bd decays yields 
B(B d -> no charm) < 14% at 2a [13]. This would 
give B(B d -> X s t + t~) < 12.5%, by subtracting the 
SM contribution of 1.5% to charmless B decays [14]. 
The compatibility between these bounds and a large 
B(B S —> t + t~) needs to be investigated in a model- 
independent framework as well as in the context of spc- 



2 



cific models, which wc shall do in this paper. 

The models with scalar lcptoquarks [7, 10, 15] and with 
flavor-dependent Z 1 couplings [16-18] have been pro- 
posed in the literature in order to enhance the b — > st + t~ 
decay rates, to account for the dimuon anomaly as well as 
the measurements of width difference in the B s -B s sys- 
tem and the CP-violating phase B s — > J/tpcf>. We explore 
the two models above to check whether they arc still able 
to account for all the data, at the same time allowing for 
an enhanced dimuon asymmetry. 



II. MODEL-INDEPENDENT ANALYSIS: 
HAMILTONIAN FOR B s -B s MIXING 

The NP that contributes to the decay modes b — > 
ST + T - directly influences B s -B s mixing, by contribut- 
ing to the dispersive (M^ 2 ) as well as absorptive (F 12 ) 
part of the effective Hamiltonian for this mixing. These 
contributions may result in an enhanced value of the life- 
time difference Ar s and the CP-violating phase (f>^ L = 
Arg[— M^ 2 /r i2 ] and, hence, will be restricted by the re- 
cent measurements of B s — > J/ip<fi at the LHCb [19]. In 
this section, we shall explore what the data have to say 
about the NP contributions to M\ 2 and T 12 . 
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FIG. 1. The goodness-of-fit contours in the IM^Hr^fl 
plane, where the phases of Mi-f and Fi2 are varied over. 
The red, blue, green and brown regions have x values less 
than 2.3, 6.18, 11.83, and 19.35, respectively, corresponding 
to regions allowed at la, 2a, 3a, and 4<r, respectively. The 
same convention is followed in all figures. 



Wc follow the recent analysis [20], with the inclusion 
of newly available data. We parameterize the NP by 



-pSM , -pp 
1 12 T 1 1 



NP 



(4) 



to 



0.04i 



where we take the central values of M 12 and T' 
be Mff = 8.65 c^" - 04 )' ps" 1 and Tf^ 1 = 0.0435 e" 
ps _1 . We investigate constraints on the two complex 
quantities (four parameters) M^ 2 P and from the 

measured values of Am s , AT S , and the CP-violating 



phase in the B s — > J/ip4> decay 



). The results of 

the x 2 fit projected in the (IM^I-Ir^)) plane arc shown 



s 



in Fig. 1. Wc assume all the measurements to be indepen- 
dent for simplicity, add the theoretical and experimental 
errors in quadrature, and vary over the phases of the two 
quantities Mf p and rf 2 p . 

The origin in this figure is the SM, which has \ 2 — 15. 
This dramatically quantifies the failure of the SM to ac- 
commodate all the current data. There are two regions in 
the parameter space that are consistent with the current 
data, one of which is also consistent with = 0, i.e. 
there is no need for NP to contribute to the dispersive 
part of the Hamiltonian. On the other hand, consistency 
with the data even to 3cr (i.e. \ 2 < 11.83) seems to 
require a nonzero NP contribution r^ 2 p . 

The 2a preferred range of Ir^l is (0.05,0.25) ps" 1 . 
Even the lower end of this range yields B{B S — > t + t~) w 
15% if r^ 2 p is coming entirely from B s — > t + t~. Thus, 
the model-independent constraints from B s -B s mixing on 
B(B S —> t + t~) are rather weak and easily allow values 
as large as ~ 15%. 



III. MODEL INDEPENDENT ANALYSIS: 
HAMILTONIAN FOR b sr+r DECAY 

Within the SM, the effective Hamiltonian for the 
quark-level transition b — > st + t~ is 

4G £ 



t/SM 
rt off 



V2 



V*V tb x 



(5) 



D 

Ci(n)Oi(n) m b {sa^P R b) 



9 ^7 (sYP L b) f^r + C w ^ (s^P L b) f 7Al75 r }, 



+ C, 



where Pl.r = (1 =F 7s)/2- The operators d (i = 1, ..6) 
correspond to the Pi of Ref. [21], and mj = to&(/x) is the 
running &-quark mass in the MS scheme. We use the SM 
Wilson coefficients as given in Ref. [22]. 

We now parameterize NP through the addition of new 
operators with distinct Lorentz structures to the effective 
Hamiltonian for b — > st + t~. Wc consider only scalar- 
pscudoscalar and vector-axial vector operators and ex- 
clude tensor ones. The new effective Hamiltonian is 

H cS (b -> sT+r-) = Hlf + + H C SP , (6) 



where the new operators are 

2 

?Cff = - 9 -^[Cv (s^PLb) + C'y {s^P R b)\ (t 7m t) 

-^§-[C A (irfP L b) + C' A (s^P R b)] (t 7m 75t) 
Ktl = -%-[C S (sP R b) + C' s (sP L b)] (fr) 



A 2 

g 2 

A 2 



[C P (sP R b) + C'p (sP L b)} (f 75 r) 



(7) 



In the above expressions, g NP and A are the NP coupling 
constant and mass scale, respectively The Ci,G[ (i — 
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V, A, S, P) are the unknown NP Wilson coefficients. Note 
that Cy and C' v do not contribute to the decay B s — > 
t + t~, but can contribute to the other b — > st + t~ modes. 
The indirect bounds on these NP operators from b — > sj 
and b —> s£ + £~ (I = e, fj,) were considered in [13] and 
found to be weaker than the direct bounds obtained from 
the ratio tb s /i~B d ■ 

In terms of all the Wilson coefficients, the branching 
ratio for this decay is given by 



B(B S 



G' 1 F a l m m BjBjB s L 4m? 



64tt 3 





t C s -C' s 


2 

+ 




mj, + m s 





2m T 



Cp-C' P 



[(V tb V t *)C 10 + Z(C A -C' A )] 



(8) 



where £ ee (gjjp/A )(v2/4GF)(47r/a em )- We now ex- 
amine the decays Bd — > Kt + t~, Bd — > K*t + t~ and 
Bd — > X s t + t~ in the presence of NP operators. The 
theoretical expressions are taken from [23, 24], with the 
replacement of m M with m T . We first consider one NP 
coupling at a time. We fix the value of the coupling by 
requiring that B(B S — > t+t~) ~ 3.5%. 

Note that when a single NP operator is present that 
enhances B(B S —> t + t~) to percent levels, since the SM 
contribution is negligible, the branching ratio would de- 
pend only on the magnitude of the NP coupling and 
not on its phase. The relevant branching ratios ob- 
tained by using the values for the NP couplings that yield 
B(B S -^t+t") = 3.5% are shown in Table I. 



Model 


B d -> Kt + t~ 




B d -> XsT+t- 


SM 


0.96 x 1CT 7 


1.15 x 10~ 7 


3.6 x 10~ 7 


Ca 


0.35 % 


0.19 % 


0.80 % 


C'a 


0.35 % 


0.19 % 


0.80 % 


Cs 


0.04 % 


0.01 % 


0.05 % 


C's 


0.04 % 


0.01 % 


0.05 % 


C P 


0.12 % 


0.03 % 


0.15 % 


C'p 


0.12 % 


0.03 % 


0.15 % 



TABLE I. Branching ratios of B d -> Kt + t~ , B d ->■ K*t + t~ 
and Bd — > X b t + t~ , when B(B S — > t + t~) = 3.5% and only 
one type of NP coupling is present. A kinematic cut q 2 > 
14.23 GeV 2 has been imposed in our theoretical calculations 
for the exclusive modes. We have set g NP = gr EW = 0.65 and 
A = 1 TeV. The theoretical errors on the branching ratios 
may be taken to be ~ 25%. 

From Table I we see that if there is only a single NP 
operator, the branching ratio of all these decay modes 
is less than 0.8%, and if the NP operator is scalar or 
pscudoscalar, it is even less than 0.2%. Therefore, the 
effect on the lifetime ratio tb s / TQ d will be small. In any 
case, the increase in Bd decay rate tends to increase this 



ratio and, hence, bring it closer to its central value of 
1.002. It is also clear that the current 90% upper bound 
on the BR of Bd — > Kt + t~ is consistent with B(B S — > 
t + t _ ) as long as the single NP operator is not of the 
axial vector kind. 

Now we shall consider the implications of the presence 
of more than one NP operator. This is of course a more 
realistic scenario, since in any specific NP model, more 
than one effective operator is generally produced. For 
example, if the NP is leptoquark [7, 10, 15], in general 
we expect both (S ± P) and (V ± A) operators to be 
generated. And if the NP is a flavor-dependent Z' [16], 
one obtains all the (V ± A) operators. 

The decay Bd — > Kt + t~ depends on the combina- 
tions {C A + C' A ), (C s + C's), and (C P + C' P ) [23, 24], 
while the B s — > t + t~ decay depends on the combina- 
tions [C A - C" A ), (C s - C's), and (C P - C'p). Therefore, 
it is always possible to enhance B(B S — > t + t~) with- 
out affecting Bd — > Kt + t~ significantly. On the other 
hand, Bd — > K*t + t~ depends on the same combina- 
tion of NP operators as B s — > t + t~ , except through 
the Fi term (see Eqs. (A.7) and (A. 12) of [23]) that de- 
pends on (Ca + C' A ), but always adds positively to the 
decay rate. Therefore, the increased branching ratio of 
B s — » t + t~ is expected to also increase the branch- 
ing ratio of Bd — > K*t + t~ . The branching ratio of 
Bd — > X s t + t~ will also naturally increase. Therefore, 
the latter two decay modes will be directly useful to con- 
strain B(B S — > t + t~) even when multiple NP operators 
are present. 

Since NP of the kind b — > st + t~ contributes to both 
B s and Bd decays (through B s — > t + t~ and Bd — s- 
X s t + t~ , respectively), the lifetime ratio T Ps jT Pd can be 
consistent with the observation even if B(B S — > t + t~) 
and B(Bd — > X s t + t~) are large, since their contribu- 
tions to the decay widths tend to cancel each other. The 
bound on B(Bd —> Kt + t~) also offers constraints, how- 
ever, these are only marginal, due to the arguments given 
above. Note that the consistency with the bounds on 
B s -t t + t~, B d -t X s t + t~, and B d — > Kt + t~ does 



not need fine-tuning, 
(which leads to B(B £ 



For example, for \Ca — C' A \ =3.8 



t + t~ ) = 15%), any set of values 
of \Ca + C' A \ and \C V + C v \ in the range [0,1.4] are al- 
lowed by the upper bound on B(Bd — > Kt + t~). The 
value of |CV| and \C' V \ may now be chosen to be ~ 3 
in order to make B(Bd — > X s t + t~) ~ 12%, so that the 
TB,/TB d constraint is satisfied. Note that further con- 



straints on B(Bd — > X s t + t ) 
limit on B(B S — > r + r _ ) from 



would decrease the upper 
- 15%. 



IV. LEPTOQUARK 

Leptoquarks (LQ) are particles whose quantum num- 
bers are such that they couple to both the quarks and 
the leptons of the SM. Vector leptoquarks are predicted 
in many NP models like models of grand unification based 
on SU(5) [25] and SO(10) [26, 27] while scalar lepto- 
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FIG. 2. The predictions of (<^ /V "*, Ar s ) within the scalar 
leptoquark model (Mlq = 250 GeV), overlaid on the experi- 
mental constraints through B s — > J/ipcj). 

quarks can arise in models of supersymmetry with It- 
parity violation [28-31] and in extended technicolor mod- 
els [32] where leptoquark states appear as bound states 
of tcchnifcrmions. With the pre-LHCb data, a third- 
generation scalar leptoquark was shown to provide an ex- 
planation of the dimuon anomaly [20] . In this section we 
revisit the leptoquark explanation of the dimuon asym- 
metry in light of the recent LHCb data. We consider an 
SU(2) singlet scalar leptoquark S with the Lagrangian, 

£ LQ DA ii (^).P^ + li.c. 1 (9) 

where (d c )j are the charge conjugate of the down-type 
quark fields, li are the charged lepton fields and i 7 j 
are the generation indices. The field S contributes to 
both M^ 2 P and r^ 2 p at one loop [13] and the contribution 
is proportional to the square of the effective coupling, 

n LQ = A32A33. 

A 95% C.L. lower bound of 210 GeV on the mass of 
a third-generation scalar leptoquark decaying to a b r fi- 
nal state was reported in [33]. Here, we show our results 
for the leptoquark mass Mlq = 250 GeV. It is worth 
mentioning that the above bound depends on specific as- 
sumptions and can be evaded if those assumptions are 
changed, however our conclusions do not change even if 
lower masses are considered. 

In Fig. 2 we show the prediction of the leptoquark 
model for two values of |Iilq| in the (<fii^^ , AT S ) plane, 
superposed on the recent results from the LHCb [19]. The 
phase of Iilq has been varied over. It is observed that a 
significant enhancement of AT S above the SM prediction 
is not possible in this model. 

A x 2 fit in the (Arg[liLQ], |Iilq |) parameter space, using 
the constraints on Am s [6] in addition to those on Ar s 

and <j)i'^ mentioned above, yields the maximum value 
of |h LQ | to be only around 0.05 (for M LQ = 250 GeV 
at 95% C.L.). This gives B(B S -> t+t") < 0.3%. A 
leptoquark model, thus, cannot raise to B(B S — > t + t~) 
to the level of a percent. It is, therefore, also not enough 
to explain the dimuon anomaly. In fact, we have checked 
that if Ag L is included in the fit , then the x 2 > xIm ~ 15 



even for very low values of leptoquark mass. This means 
that the LQ model cannot do better than the SM as far 
as the explanation of the dimuon anomaly is concerned. 

V. FLAVOR CHANGING Z' 

In this section we consider a flavor-changing Z' [16-18] 
gauge boson and examine whether it can offer an expla- 
nation of the dimuon anomaly while still being consistent 
with the recent LHCb data. We consider the Lagrangian 
density 

L z , D (R sb [s^P L b]Z'v + h.c.) + R TT [f 7Al P L r]Z'r (10) 

Note that Am s gets tree-level Z' contribution from Lz< 
and is proportional to (R sb ) 2 . This restricts the value of 
R sb to be very small. On the other hand, rf 2 is sensi- 
tive to the combination of couplings A = R sb R TT , and we 
have already seen in Fig. 1 that a large value of r^ 2 p is re- 
quired for a solution of the dimuon anomaly. This means 
that a large value of A will be needed in the Z' model. 
With a small R sb , this could imply a very large value of 
R TT , where the calculations would become nonperturba- 
tive and, hence, unreliable. The situation would become 
worse with increasing values of Z' mass. Therefore we 
have to stay away from the nonperturbative region. 

Lower experimental bounds on the Z' mass can be a 
serious problem. Direct searches of pair production of 
t + t~ at the Tevatron have provided a lower bound on 
the Z' mass to be 399 GeV at 95% C.L. [34]. However 
since this bound assumes SM-like couplings of Z' to all 
quarks, it can be bypassed if the Z' is assumed to couple 
very weakly (or not at all) to the first-generation quarks. 
This allows us to consider a very light Z' . with mass 
as low as M z > = 7 GeV. Note that since M z > « M Bs , 
the decay width of Z' can affect the results (at large 
couplings, the width may be as large as a few GeV) and 
has been taken into account in our calculations. 

Note that, in principle, bounds on the effective cou- 
pling in Z — > t + t~ as obtained in [35] can put correlated 
constraints on R TT and Mz> ■ However, the fit therein as- 
sumes a small imaginary part of the effective coupling, 
while the one- loop correction [36] to Z — >• t + t~ through 
a Z' exchange has a large imaginary part for small AIz> ■ 
The constraints in [35] are, therefore, not applicable in 
our case. A further analysis of Z —> t + t~ , allowing for 
an imaginary contribution to the effective coupling, may 
yield stronger constraints. 

In Fig. 3 we show the predictions of the Z' model in 
the (0s^^,Ar s ) and (^ L , AF S ) planes for some sam- 
ple values of the couplings for illustration. (Note that 
± </>sl [7 j 37].) The figure shows that there are 
values of the couplings that can be consistent with the 
J/4>4> data as well as the dimuon asymmetry data to 
within 2cr. The values of B(B S — > t + t~) in these allowed 
regions can be as high as ~ 20%. However, constraints 
from the measurements of Ts q /rs ti , B(B e i — > X s t + t~) 
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FIG. 3. The upper panel shows the predictions of the Z' 
model with sample values of the couplings in the (<pi^^ — 
AT a ) plane overlaid with the constraints from B 3 — >• J/ip(j>. 
The lower panel shows the corresponding predictions in the 
(cf>f L — AF S ) plane overlaid with the experimental constraints 
from A b SL and a d SL - We have used M z / = 7 GeV. 

and B(B c i — > Kt + t~) still continue to apply. In the case 
of a specific model such as this, the constraints would 
be more severe than those in Sec. Ill since the effective 
couplings (Cy,C' A ) are now related to each other. 

Figure 4 shows the result of a x 2 fit in the (\R sb \,R TT ) 
plane using all the experimental data: Am s , Ar s , 
Ag L ,seem consistent with all data to within 2a, the re- 
sults are reliable only when the couplings are perturba- 
tive, i.e. a RTT = (i? TT ) 2 /47r < 1. We indicate this value 
of R TT by a horizontal line in the figure, and consider 
only the region below this line as the valid one. (The 
region below the line but very close to it may still have 
significant higher-order corrections.) 

Even by demanding pcrturbative couplings in addi- 
tion to x 2 < 6.18 (i.e. 2a level), the branching ra- 
tio B(B S —> t + t~) can be as high as ~ 20%. Thus, 
the Am s , J / i\)<\), and A^ L data by itself cannot put 
strong constraints on the value of this branching ra- 
tio. The strongest constraints here come from tb s j T~g d 
and B(Bd — > Kt + t~). Taking into account these con- 
straints, the value of B(B S — > t + t~) can be as high as 
5%. This can happen, for example, with R sb = 15 x 1CP 6 , 
R TT = 3.0, which makes B(B d ~> X s t+t-) tn 1.5% and 




20 30 

|R sb | X 10 6 

FIG. 4. Goodness-of-fit contours in the (\R sb \-R TT ) plane for 
M z / = 7 GeV. The horizontal line corresponds to R TT = y/4ir. 
Contours corresponding to B(B S — > t + t~) values of 3.5%, 
7%, 12%, and 20% are also shown. 

B(Bd — > Kt + t~) rj 0.35%, the former allowing a relax- 
ation of the constraint from tb s /tb^ ■ 

Note that this value of B(B S — > t + t~) can alleviate 
the dimuon anomaly (from \ 2 = 15.1 for the SM, to 
X 2 = 9.5); however, it is not enough to explain it. Some 
contribution from NP in the Bd — Bd mixing may also be 
needed, as was recently conjectured in [38]. 



VI. CONCLUDING REMARKS 

We have investigated the possible enhancement of 
B(B S — > t + t~) that may help explain the observed 
anomalous like-sign dimuon asymmetry Taking into ac- 
count the constraints from the lifetime ratio of B s and 
Bd, as well as the measurements of related b — > st + t _ 
decay modes, we find that an enhancement up to B(B S — > 
t + t~) ~ 15% is allowed at 2c. This bound may decrease 
with further constraints on B{Bd X s t + t~). 

Within the context of specific models, an enhancement 
of B(B S — > t + t~) to more than 0.3% is not possible with 
the leptoquark model, but the model with a light flavor- 
changing Z' that does not couple to light quarks can 
increase it to 5%. This helps alleviate the Ag L anomaly 
to some extent, but cannot account for it entirely, and 
contribution from NP in the Bd sector may be needed. 

Our results with leptoquarks are similar to those in 
[13], however the updated LHCb data used by us has 
made the constraints on leptoquark parameters even 
stronger. On the other hand, the light Z' employed by us 
allows much larger values of the branching ratio than that 
predicted therein. Accounting for the dimuon anomaly to 
within la, as indicated in [16], is not possible even with 
a light Z' , because of the new stronger LHCb constraints 
and the requirement of pcrturbative couplings. 

If B(B S — > t+t~) is at the percent level, it will soon 
be within the reach of experiments, and could play an 
important role in our understanding of physics beyond 
the SM. The direct measurement of this branching ratio 
should, hence, be a high-priority. 



6 



ACKNOWLEDGMENTS S. Nandi, and S. K. Patra for valuable discussions and 

comments. 

We thank A. K. Alok, C. Bobeth, R. Godbole, M. 
Grunewald, U. Haisch, A. Kundu, A. Lcnz, D. London, 



[1] V. M. Abazov et al. [DO Collaboration], Phys. Rev. D 
82, 032001 (2010) [arXiv: 1005.2757 [hep-ex]]; 

[2] V. M. Abazov et al. [DO Collaboration], Phys. Rev. Lett. 
105, 081801 (2010) [arXiv: 1007.0395 [hep-ex]]. 

[3] V. M. Abazov et al. [DO Collaboration], Phys. Rev. D 84, 
052007 (2011) [arXiv:1106.6308 [hep-ex]]. 

[4] CDF Collaboration, CDF Note 9015 (2007) 

http : //www-cdf . f nal . gov/physics/new/bottom/ 
070816 .blessed-acp-bsemil 

[5] A. Lenz, U. Nierste, J. Charles, S. Descotes-Genon, 
A. Jantsch, C. Kaufhold, H. Lacker, S. Monteil, V. Niess 
and S. T'Jampens, Phys. Rev. D 83, 036004 (2011) 
arXiv: 1008. 1593 [hep-ph]; A. Lenz and U. Nierste, 
arXiv:1102.4274 [hep-ph]; A. Lenz, Proceedings of the 
CKM workshop 2010 in Warwick, arXiv:1205.1444 [hep- 
ph]. 

[6] D. Asner et al. [Heavy Flavor Averaging Group], 

arXiv:1010.1589 [hep-ex], and updates online at 

http : //www. slac . Stanford. edu/xorg/hf ag. 
[7] A. Dighe, A. Kundu and S. Nandi, Phys. Rev. D 82, 

031502 (2010) [arXiv:1005.4051 [hep-ph]]. 
[8] C. W. Bauer and N. D. Dunn, Phys. Lett. B 696, 362 

(2011) [arXiv:1006.1629 [hep-ph]]. 
[9] R. Aaij et al. [LHCb Collaboration], Phys. Rev. Lett. 

108, 101803 (2012) [arXiv:1112.3183 [hep-ex]]. 
[10] Y. Grossman, Z. Ligeti and E. Nardi, Phys. Rev. D 55, 

2768 (1997) [hep-ph/9607473] . 
[11] K. Flood (BaBar Collaboration), Proc. Sci. 

ICHEP2010, 234 (2010). 
[12] D. Buskulic et al. [ALEPH Collaboration], Phys. Lett. B 

343, 444 (1995). 
[13] C. Bobeth and U. Haisch, arXiv: 1109. 1826 [hep-ph]. 
[14] A. Lenz, U. Nierste and G. Ostermaier, Phys. Rev. D 56, 

7228 (1997) [hep-ph/9706501]. 
[15] S. Davidson, D. C. Bailey and B. A. Campbell, Z. Phys. 

C 61, 613 (1994) [hep-ph/9309310]. 
[16] A. K. Alok, S. Baek and D. London, JHEP 1107, 111 

(2011) [arXiv:1010.1333 [hep-ph]]. 
[17] X. -Q. Li, Y. -M. Li, G. -R. Lu and F. Su, JHEP 1205, 

049 (2012) [arXiv:1204.5250 [hep-ph]]. 
[18] H. D. Kim, S. -G. Kim and S. Shin, arXiv: 1205.6481 [hep- 
ph]. 



[19] LHCb Collaboration, LHCb-CONF-2012-002 (2012) 
[20] A. Dighe, D. Ghosh, A. Kundu and S. K. Patra, Phys. 

Rev. D 84, 056008 (2011) [arXiv:1105.0970 [hep-ph]]. 
[21] C. Bobeth, M. Misiak and J. Urban, Nucl. Phys. B 574, 

291 (2000) [arXiv:hep-ph/9910220]. 
[22] S. Descotes-Genon, D. Ghosh, J. Matias, M. Ramon, 

JHEP 1106, 099 (2011). [arXiv:1104.3342 [hep-ph]]; 

Proc. Sci. EPS -HEP2011, 170 (2011) [arXiv: 1202.2172 

[hep-ph]]. 

[23] A.K.Alok et.al, JHEP 1002, 053 (2010) [arXiv:0912. 1382 
[hep-ph]]; 

[24] A.K.Alok et.al, JHEP 1111, 121 (2011)[arXiv:1008.2367 
[hep-ph]]; JHEP 1111, 122 (2011) [arXiv:1103.5344 [hep- 
ph]]. 

[25] H. Georgi and S. L. Glashow, Phys. Rev. Lett. 32, 438 

(1974) . 

[26] H. Georgi, AIP Conf. Proc. 23, 575 (1975). 

[27] H. Fritzsch and P. Minkowski, Annals Phys. 93, 193 

(1975) . 

[28] G. R. Farrar and P. Fayet, Phys. Lett. B 76, 575 (1978). 
[29] S. Weinberg, Phys. Rev. D 26, 287 (1982). 
[30] N. Sakai and T. Yanagida, Nucl. Phys. B 197, 533 (1982). 
[31] C. S. Aulakh and R. N. Mohapatra, Phys. Lett. B 119, 
136 (1982). 

[32] E. Farhi and L. Susskind, Phys. Rept. 74, 277 (1981). 
[33] V. M. Abazov et al. [DO Collaboration], Phys. Rev. Lett. 

101, 241802 (2008) [arXiv:0806.3527 [hep-ex]]. 
[34] D. Acosta et al. (CDF Collaboration), Phys. Rev. Lett. 

95, 131801 (2005) [hep-ex/0506034]. 
[35] ALEPH, DELPHI, L3, OPAL, SLD and LEP Elec- 

troweak Working Group and, SLD Electroweak Group 

and SLD Heavy Flavour Group Collaborations, Phys. 

Rept. 427, 257 (2006) [hep-ex/0509008] . 
[36] U. Haisch and S. Westhoff, JHEP 1108, 088 (2011) 

[arXiv:1106.0529 [hep-ph]]. 
[37] A. Lcnz, Nucl. Phys. Proc. Suppl. 177-178, 81 (2008) 

[arXiv:0705.3802 [hep-ph]]. 
[38] A. Lenz, U. Nierste, J. Charles, S. Descotes-Genon, 

H. Lacker, S. Monteil, V. Niess and S. T'Jampens, Phys. 

Rev. D 86, 033008 (2012) arXiv: 1203. 0238 [hep-ph]. 



